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Government Funded Invention 
The Government may have certain rights to the invention as provided by the terms of 
Agreement DTFAO 1-98-00031 awarded by the Federal Aviation Administration. 

Background of the Invention 

7. Field of the Invention 

The invention is related to the field of radar systems, and in particular, to radar systems 
that transfer pulses at a target to determine the range and velocity of the target. 

2. Statement of the Problem 

A Doppler weather radar system determines the range and velocity of atmospheric 
structures, such as wind and precipitation. The Doppler weather radar system transfers pulse 
trains toward the target structure and receives echoes reflected off of the target structure. The 
Doppler weather radar system processes the echoes to determine the range and velocity of the 
atmospheric structure. In addition, the Doppler weather radar system determines other variables, 
such as differential reflectivity, differential phase, linear depolarization ratio, and co-polar 
correlation coefficient. 

To determine the target range, the radar system determines the time that it takes for a 
pulse to propagate to the target and for its echo to propagate back to the radar system. The range 
is then determined based on the expected speed of the pulses/echoes and the elapsed time. 
However, the echoes from consecutive pulses may be confused by the radar system. For 
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example, a first pulse may reflect off of the target to provide a first echo, and a second pulse may 
reflect off of the target to provide a second echo. However, the first pulse may continue to 
propagate and reflect off of another structure that is further away from the radar system than the 
target structure to produce a third echo. The third echo may arrive at the receiver around the 
5 same time as the second echo. The receiver may confuse the two echoes and provide an 

ambiguous or inaccurate target range. These confounded echoes are referred to as range overlaid 
echoes in the literature. 

To counter the ambiguous range problem, it is desirable to separate consecutive pulses by 
a relatively long Pulse Repetition Time (PRT). The longer PRT allows all echoes from the first 
10 pulse to clear the radar system before the echoes from the second pulse arrive. This prevents the 
receiver from confusing the echoes of consecutive pulses. 

To determine target velocity, the radar system determines the phase shift that occurs 
between two echoes returned from a common target, since the velocity of the target affects this 
phase shift according to the Doppler effect. To get proper phase shift data, the target must be 
1 5 sampled at the Nyquist rate or at a higher rate, and these rates require a shorter PRT to 

adequately sample for unambiguous target velocity. 

The velocity determination can be further complicated by the introduction of phase shift 
that is not attributed to the velocity of the target structure. For example, additional phase shift is 
introduced to both the pulse and its echo by the particles that are in between the radar system and 
20 the target structure. Since this other phase shift is not the Doppler phase shift of interest, the 
other phase shift introduces ambiguity into the velocity determination. 

Unfortunately, larger unambiguous range determination requires a longer PRT, but larger 
unambiguous velocity determination requires a shorter PRT. This problem of simultaneously 
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obtaining larger unambiguous range and larger unambiguous velocity is termed in the literature 
as the Doppler Dilemma. It is expressed by the following equation: 
r a v a = cA/8; 

where r a is the unambiguous range, v a is the unambiguous velocity, c is the speed of light, and X 
is the wavelength of the transmitted radar signal. Various techniques have been tried to 
simultaneously extend both the unambiguous range and velocity of the target, despite the 
conflicting pulse repetition requirements. 

One technique to handle the range/velocity ambiguity problem is to use dual scans where 
the first scan has a short PRT for velocity determination, and the second scan has a longer PRT 
for range determination. Unfortunately, the dual scanning technique may take an undesirable 
amount of time to perform both scans, and range overlaid echoes can contaminate the short PRT 
scan. 

Another technique to handle the range/velocity ambiguity problem is to alternate the 
period between the pulses between a short PRT and a long PRT, called a staggered PRT. The 
unambiguous velocity can be extended to an equivalent uniform PRT of the difference of the 
long and short PRT periods, while the unambiguous range can be extended to an equivalent 
uniform PRT of the sum of the long and short PRT. Unfortunately, complex filters are required 
to remove noise caused by ground clutter when the PRT is alternated. 

Another technique to mitigate range-velocity ambiguity is to use systematic phase coding 
of the transmitted pulses. In this way, the various range overlaid echoes can be separated using 
the phase code information. This technique works well for two overlaid echoes, but becomes 
very difficult when three or more echoes are overlaid. 
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Another technique to handle the range/velocity ambiguity problem is to transfer a first 
pulse train at a first polarity and with a set pulse repetition time, and to transfer a second pulse 
train at a second polarity, but with an alternating PRT. The alternating PRT means that, for the 
second pulse train, a first pulse is transmitted slightly before a corresponding first pulse in the 
5 first pulse train, and a second pulse is transmitted slightly after a corresponding second pulse in 

the first pulse train. Thus, pulses in the second pulse train alternate between transmission 
slightly before or after pulses in the first pulse train. Unfortunately, the continually alternating 
PRT adds undesired complexity to ground clutter filtering and to phase coding. 

10 Summary of the Solution 

Examples of the invention include a radar system and its method of operation. The radar 
system comprises a transmission system and a reception system. The transmission system is 
configured to transfer a first series of pulses and a second series of pulses. The first series of 
pulses and the second series of pulses have orthogonal polarizations. The first series of pulses 

15 and the second series of pulses have a same pulse repetition time. The first series of pulses and 

the second series of pulses are offset by a time amount. A target reflects energy from the first 
series of pulses to generate a first series of echoes. The target reflects energy from the second 
series of pulses to generate a second series of echoes. The reception system is configured to 
process the first series of echoes and the second series of echoes to determine a range and a 

20 velocity of the target. 

In some example of the invention, the first series of pulses and the second series of pulses 
each comprise blocks of the pulses, and the transmission system is configured to alternate the 
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blocks from the first series of pulses between a lead position and a lag position with respect to 
the blocks from the second series of pulses. 

In some example of the invention, the transmission system is configured to phase code 
the first series of pulses and the second series of pulses, and the reception system is configured to 
decode the phase coding to distinguish ones in the first series of echoes from adjacent ones in the 
second series of echoes. 

In some example of the invention, the orthogonal polarizations are substantially parallel 
and perpendicular to earth's surface. 

In some example of the invention, the reception system is configured to co-process ones 
in the first series of echoes with adjacent ones in the second series of echoes to determine the 
velocity of the target. 

In some example of the invention, the reception system is configured to co-process 
adjacent ones in the first series of echoes and co-process adjacent ones in the second series of 
echoes to determine the range of the target. 

In some example of the invention, the pulse repetition time is in a range from one-half 
millisecond to four milliseconds. 

In some example of the invention, the time amount is in a range from one microsecond to 
300 microseconds. 

In some example of the invention, the target comprises an atmospheric structure. 

In some example of the invention, the radar system comprises a Doppler weather radar 

system. 
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Description of the Drawings 
The same reference number represents the same element on all drawings. 
FIG. 1 illustrates a radar system in an example of the invention. 
FIG. 2 illustrates a first and second series of pulses in an example of the invention. 
FIG. 3 illustrates a first and second series of pulses in an example of the invention. 
FIG. 4 illustrates radar system processing in an example of the invention. 

Detailed Description of the Invention 

FIGS. 1-4 and the following description depict a specific example to teach those skilled 
in the art how to make and use the best mode of the invention. For the purpose of teaching 
inventive principles, some conventional aspects have been simplified or omitted. Those skilled 
in the art will appreciate variations from this example that fall within the scope of the invention. 
Those skilled in the art will appreciate that the features described below can be combined in 
various ways to form multiple variations of the invention. As a result, the invention is not 
limited to the specific example described below, but only by the claims and their equivalents. 

FIG. 1 illustrates radar system 100 in an example of the invention. Radar system 100 
determines the range and velocity of target 105. Target 105 could be an atmospheric structure, 
such as wind or precipitation. Radar system 1 00 could be a Doppler weather radar system, Lidar 
system, or other radar-like system. 

Radar system 100 includes transmission system 101 and reception system 102. 
Transmission system 101 receives control signal 103. Transmission system 101 transfers a first 
series of pulses 1 1 1 and a second series of pulses 121 toward target 105. Target 105 reflects 
energy from first series of pulses 1 1 1 to generate first series of echoes 112. Target 105 reflects 
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energy from second series of pulses 121 to generate second series of echoes 122. Reception 
system 102 receives and processes first series of echoes 112 and second series of echoes 122 to 
determine the range and velocity of target 105. Reception system 102 provides output signal 104 
that indicates the range and velocity of target 1 05. 

Other variables may also be determined by processing the echoes, such as the differential 
reflectivity, differential phase, linear de-polarization ratio, and co-polar correlation coefficient. 
The differential reflectivity is the co-polar power differential between echoes 112 and 122. The 
differential phase is the co-polar phase differential between echoes 1 12 and 122. The linear de- 
polarization ratio is a ratio of the cross-polar power to the co-polar power. The co-polar 
correlation coefficient indicates the degree of similarity between echoes 112 and 122. 

FIG. 2 illustrates first series of pulses 1 1 1 and second series of pulses 121 in an example 
of the invention. On FIG. 2, the horizontal axis represents time and the vertical axis represents 
magnitude and polarity. First series of pulses 1 1 1 includes pulses 211-214. Second series of 
pulses 121 includes pulses 221-224. First series of pulses 1 1 1 and second series of pulses 121 
both have the same pulse repetition time Tp. First series of pulses 111 and second series of 
pulses 121 are offset from one another by an offset time amount To , where the first series of 
pulses 1 1 1 leads the second series of pulses 121 by the offset time amount To. 

In some examples, the pulse repetition time Tp is in a range from one-half millisecond to 
four milliseconds, and the time amount To is in a range from one microsecond to 300 
microseconds. In addition, control signal 103 may specify new values for To and T P , and in 
response, transmission system 101 would change the prior values of T 0 and T P to the newly 
specified values for the transmission of pulse series 1 1 1 and 121. 
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First series of pulses 1 1 1 and second series of pulses 121 have orthogonal polarizations. 
On FIG. 2, the pulses above the time axis have a first polarization, and the pulses below the time 
axis have a second polarization - where the first polarization is orthogonal to the second 
polarization. For example, the polarization of pulse 21 1 is orthogonal to the polarization of pulse 
221. In some examples, one of the orthogonal polarizations is substantially perpendicular to the 
earth's surface, and the other one of the orthogonal polarizations is substantially parallel to the 
earth's surface. These perpendicular and parallel polarizations are respectively referred to as 
horizontal and vertical polarizations. 

To determine the range of target 105, reception system 102 co-processes adjacent echoes 
from first series of echoes 112 and co-processes adjacent echoes from second series of echoes 
122. For example, reception system 102 would co-process echoes for adjacent pulses 21 1 and 
212 in first series of echoes 111, and co-process echoes for adjacent pulses 221 and 222 in 
second series of echoes 122. The term "co-process" means that data from the echoes are used in 
the same equations at the same time. 

Advantageously, the relatively long time Tp between the adjacent pulses of the same 
series provides a relatively slow sampling rate to unambiguously determine the target range. The 
slow sampling rate allows echoes from the prior pulse to clear before echoes from the next pulse 
in the same series are received. Thus, reception system 102 is unlikely to confuse echoes from 
adjacent pulses in the same series. 

To determine the velocity of target 105, reception system 102 co-processes adjacent 
echoes from first series of echoes 112 and second series of echoes 122. For example, reception 
system 102 would co-process echoes for adjacent pulses 211 and 221 and co-process echoes for 
adjacent pulses 212 and 222. 
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Advantageously, the relatively short time amount To between the adjacent pulses of the 
different series provides a relatively fast sampling rate to unambiguously determine the target 
velocity. The fast sampling rate can handle the Doppler frequency shifts of high velocity targets. 

Reception system 102 must also distinguish between echoes from the adjacent pulses of 
the different series where the separation between the adjacent pulses is the short time amount To. 
For example, echoes from pulses 211 and 221 must be distinguished. The different polarities of 
the pulses from the different series can be used to distinguish their respective echoes. To handle 
cross-polar interference between these pulses, transmission system 101 phase codes first series of 
pulses 1 1 1 differently from second series of pulses 121 . Reception system 102 decodes this 
phase coding to distinguish first series of echoes 112 from second series of echoes 122. For 
example, reception system 102 uses the phase coding to distinguish the echoes from adjacent 
pulses 21 1 and 221, or to distinguish the echoes from adjacent pulses 212 and 222. 

FIG. 3 illustrates first series of pulses 1 1 1 and second series of pulses 121 in an example 
of the invention. FIG. 3 is a continuation of FIG. 2, starting with pulses 214 and 224 that appear 
on both FIGS. 2-3. The second series of pulses 121 continues to use the standard pulse repetition 
time T P , so pulse 225 follows pulse 224 by T P , and pulse 226 follows pulse 225 by T P . However 
in the first series of pulses 1 1 1, pulse 215 does not follow pulse 214 by the standard pulse 
repetition time Tp, but instead, is delayed by an additional 2To. Starting with pulse 216 (which 
follows delayed pulse 215), the first series of pulses 1 1 1 again uses the standard pulse repetition 
time T P . Note that until pulses 214 and 224, the first series of pulses 111 led the second series of 
pulses 121 by the time offset To, but after pulses 214 and 224, the first series of pulses 1 1 1 lags 
the second series of pulses 121 by the time offset To. After a set number of pulse pairs (two 
pulses of different series that are offset by the time amount To), the above-described pulse delay 
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process will occur for the second series of pulses 121, so that the first series of pulses 1 1 1 will 
again lead the second series of pulses by the offset time To. This process repeats, so that the first 
series of pulses 1 1 1 and the second series of pulses alternate between leading and lagging one 
another by the offset time To. 

The lead/lag positioning alternates after a number of pulse pairs, where the number could 
be in a range of 16-256 or 32-128. Control signal 103 may specify the number of pulse pairs that 
occur between lead/lag alternation, and in response, transmission system 101 would change the 
prior number of pulse pairs to the newly specified number of pulse pairs for the transmission of 
pulse series 111 and 121. 

The first series and second series of pulses can each be viewed as blocks of consecutive 
pulses, such as blocks of 64 pulses. For the first block from each series, the block for the first 
series leads the block for the second series by the offset time To. For the second block from each 
series, the block for the second series leads the block for the first series by the offset time To. For 
the third block from each series, the block for the first series leads the block for the second series 
by the offset time To. This process repeats, so that consecutive blocks from the same series 
alternate between the lead and lag position with respect to the blocks of the other series. 

Advantageously, alternating the lead/lag positioning of the two series of pulses allows 
reception system 102 to separate phase shift attributed to the velocity of the target structure (the 
Doppler phase shift) from other phase shift, such as the other phase shift that is introduced by 
particles in between the radar system and the target structure. Alternating the lead/lag 
positioning changes the sign of the Doppler phase shift while leaving the sign of the other phase 
shift alone, and thus, the Doppler phase shift may be mathematically separated from the other 
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phase shift. The ability to separate the Doppler phase shift from the other phase shift improves 
the accuracy of the velocity determination. 

Even though the lead position alternates between the two pulse series, note that the pulse 
repetition time T P for the pulses in both series largely remains the same. Only a single pulse is 
periodically delayed to alternate the lead position between the two series. Advantageously, the 
consistent pulse repetition time simplifies ground clutter filtering and phase coding. In contrast, 
ground clutter filtering and phase coding are relatively complex if the pulse repetition time is 
allowed to alternate on a pulse-by-pulse basis as is done in some prior systems. 

FIG. 4 and the following section describe the processing that is used by Radar system 100 
in one example of the invention to determine the range and velocity of the target, in addition to 
other variables. Note that other processing variations could be used in other examples of the 
invention. 
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1. Staggered H and V Pulse Trains 

1 . 1 Estimation of copolar and crosspolar signals via phase coding 

First the copolar and crosspolar components of the received horizontal and vertical 
voltages are estimated via phase coding. Refer to Figure 4. At time t s the received voltages are: 



where H and V are the received voltages for the horizontal and vertical polarizations, 
respectively, and where the letter pairs refer to the received and transmitted polarizations. For 
example, vh is transmit horizontal polarization receive vertical polarization. 0, and are the 
phase codes applied to the horizontal and vertical transmitted pulses respectively. The letter i 
represents the time series index for the radar copolar resolution volume defined at time t s? and j is 
the square root of -1 . The superscript (1) refers to the samples obtained when the V pulses lag 
the horizontal pulses by tj. The superscript (2) refers to the samples obtained when the vertical 
pulses lead the horizontal pulses by t\. In the equations below, it is assumed that ti < U < T - ti; 
however, the processing principles are identical. 

Before processing the samples to retrieve the desired radar measurables, the copolar and 
cross polar time series first need to be estimated. Note that in prior systems (Sachidananda and 
Zrnic), the objective was to estimate only the radar moments of power, velocity and spectrum 
width. Here we need to obtain the copolar time series so that differential phase, velocity and 



= vv^ii^expij^ + vh^ii.ts + t^expijcPi} 
HW(t) = W 2 '(z,t 5 )exp{# 2 } + hv^^ts+tJexpijiPi} 



(1) 

(2) 
(3) 
(4) 
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copolar correlation coefficient can be estimated. The phase codes should be chosen so that 
overlaid echoes have a unique phase code. For example, if the sequence length is 64 and the 
SZ(8/64) phase code is chosen, then a good choice is = 01-3. 

It is possible that second trip echoes may overlay the received voltages given in equations 
1, 2 5 3 and 4. Under certain conditions, it is possible to separate three or more overlaid signals 
but in the following discussion it is assumed that only first trip echoes are present, i.e., all 
significant targets lie within the range cT/2 where c is the velocity of light and Tis the pulse 
repetition time. 

Recovery of the weaker echo for overlaid SZ phase coded signals consists, of the 
following steps: 

1. Cohere to the stronger signal 

2. Apply an aggressive window function to the resulting time series such as the Harming 
window. 

3. Apply a Fast Fourier Transform (FFT) to obtain the frequency spectrum. 

4. Notch out the strong signal spectrum. 

5. Inverse Fast Fourier Transform (IFFT) resulting spectrum. 

6. Cohere to the weak signal 

The time series for the strong and weak signal components of the horizontal received 
signal are designated Sh and Wh, respectively, while the time series for the strong and weak 
signal components of the vertical received signal are designated Sv and Wy, respectively. These 
time series may be either copolar or crosspolar. 
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The strong signal is estimated as: 

S%\i.,t s ) = HW(i,t s ) (5) 
S$\i,t.) = W(i,t s ) (6) 



or 



Sg>(t,t.) = HW(i,t.) • (7) 

: SP(i,t s ) = V< 2 >(t,t.) (8) 

The weak signal time series is recovered via the phase coding recovery algorithm 
outlined above. After the above processing, time series representing the copolar and crosspolar 
received signal for both the horizontal and vertical receive channel are available at all times, t s . 

The recovered horizontal copolar and crosspolar time series are now designated as, 
H (1) (i, t s ),1lV (l) (i, t s ), respectively, and the recovered vertical copolar and crosspolar time series 
are now designated as, (i, t s ), respectively. For the second block of pulses, the 

superscript (2) replaces the superscript (1). 

1 2 Calculating the radar variables 

The processing steps for the proposed pulsing scheme are: 



1. Separate the copolar and cross polar time series, as outlined above. 
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2. Calculate the copolar powers Ph, P v and crosspolar powers, P V h, Phv for time t s as: 



1=1 



Pv = _Lf vu + v® (io) 



= (12) 



3. Correct the strong signal power by subtracting the weak signal power if 10 log P s /P w < 
20 dB where P S5 P w are the strong and weak signal powers. 

4. Calculate (copolar differential reflectivity) and LDR h (linear depolarization ratio 
for horizontal transmit polarization) as: 



^ dr ~ 7T (13) 



LDR = — 
Ph 



5. Calculate the copolar differential phase and velocity phase: 



(14) 



i=l 
N 



02 = d>d P + <t> ve i = arctan^f// 2) ^;^ (16) 

i—1 



<Pvel = (02-9l)/2 



(17) 



= (6 2 + 0i)/2 (18) 
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where * represents complex conjugation. 



6. If t] is short, the standard deviation of this above velocity phase estimate, 0 ve i, can 
be high, but the associated Nyquist folding velocity will be high as compared to the 
Nyquist folding velocity associated with the time interval T. Thus, velocity phase 
is also calculated using the first lag autocovariance of the copolar vertical polarized 
samples as, 



(pveijr ~ arctan < 



E V5 (1) vff? 



k L i=l 



+ 



N-l 

E 



T V(2) T ^*(2) 
x v i+l 



(19) 



This velocity phase 0 ve i_T will have a lower Nyquist folding velocity than the velocity 
phase estimate 0 ve i but the standard deviation of the velocity phase estimate 0 ve i_T will 
in general be lower than the standard deviation of the velocity phase estimate. 0 ve i for 
small ti. Therefore, the velocity phase estimate 0 ve i can be used to place the velocity 
phase estimate 0 ve i_T in the correct Nyquist interval thereby by giving a velocity phase 
estimate with a high Nyquist phase and a low standard deviation. 



7. To calculate the copolar correlation coefficient, magnitude deconvolution should be applied to 
the recovered weak signal time series that resulted after step 6 of the phase coding outline in 
section 1.1 above. After applying magnitude deconvolution to the copolar time series, if needed, 
the copolar correlation coefficient is calculated as, 



Phv — 



Ef =1 vPh*™ 



/2 (20) 
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If the time U is small so that there is very small decorrelation between the H\ and V x 
signals, the above estimate of ph v in equation 20 can consider the zero lag correlation coefficient 
5 between the H ; and Vj time series with very little error. 

CLAIMS: 
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